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Abstract. We investigate an initial-boundary value problem for the quasi- 
linear Westervelt equation which models the propagation of sound in 
fluidic media. We prove that, if the initial data are sufficiently small 
and regular, then there exists a unique global solution with optimal Lp- 
regularity. We show furthermore that the solution converges to zero at 
an exponential rate as time tends to infinity. Our techniques are based 
on maximal Lp-regularity for abstract quasilinear parabolic equations. 
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1. Introduction and Notations 

The aim of this paper is to enhance the mathematical understanding of the 
Westervelt equation, which was proposed and analyzed in niiniin]: 

u'^^ - c^A^u' - bA^u't = k{u'%f ( 1 - 1 ) 

Here u'(t,x) = u{t,x) — uq denotes the acoustic pressure fluctuation at 
time t and position x. Furthermore, c > 0 denotes the velocity of sound, 
5 > 0 the diffusivity of sound and k > 0 the parameter of nonlinearity. 
This equation is used to describe the propagation of sound in fluid media. It 
can be derived from the balances of mass and momentum (the compressible 
Navier-Stokes equations for Newtonian fluids) and a state equation for the 
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pressure-dependent density of the fluid. A generalization of (HH) is given by 
Kuznetsov’s equation: 

u'tt - - bA^u't = k{u‘^)tt + |Vu'|(j. (1.2) 

Here the velocity fluctuation v'{t,x) = v{t,x) — vq is related to the pressure 
fluctuation by means of an acoustic potential 'ijj{t,x), such that u' = poipt-, 
v' = —V^/:. This equation is used as a basic equation in nonlinear acoustics, 
see [a nans]- We refer to [12] for a derivation of Kuznetsov’s equation. 

To avoid any confusion, we mention that in the engineering literature 
laiiaiii] also the following equation is termed Westervelt equation. 

'^tt ~ c ^xU H 2^ttt ~ 

This equation is derived in (a, for instance. However, we may still view (ED 
as a simplification of ED- 

Let H C M" be a bounded C^-domain, let J = (0,T), T S (0,oo] and 
let uq, ui be given functions on H. Then we consider the following initial¬ 
boundary value problem for the Westervelt equation. 

{ Utt — c^Au — bAut = k{u^)tt in J x H, 

= 0 in J X dil (1.3) 

{u{0),ut{0)) = (uo,ui) in H. 

Here u : J x H R is the unknown function and {uo,ui) : H —j- R^ 
are the given initial data. In this paper we prove existence and unique¬ 
ness of global strong solutions for small initial data in anisotropic Sobolev- 
Slobodeckii spaces: 

Theorem 1.1. Let n € N, p > max{n/2, n/4-|- 1}, p ^ 3/2 and suppose that 
the initial data (uo,ui) satisfy the regularity and compatibility conditions 

uo€ {v e Wp{D,) : -clan = 0}, 

^ ^ f Wp'-"/P(H) , tfp < 3/2, (1.4) 

\{ue Wp-"/^(H) :u|an = 0} , */p > 3/2. 

Then for every T € (0, cx)] there is p > 0 such that the smallness condition 

II^OIIW^ + ||Ul||^,^;-2-2/p < P 
implies that problem (ED admits a unique solution 

u G W2(0,T;Lp(H)) nM/)(0,r; W2(H)). (1.5) 

The solution map (uo,ui) u is a local isomorphism in these spaces. 

This is an extension of Theorems 1.1 and 1.2 of m , where only the case 
p = 2, n G {1,2,3} and ||mo||w 2 + l|wi|lw| < P is considered. Our result also 
implies that the regularity class in m Theorem 1.1] is not optimal. In order 
to provide optimal regularity results for the solutions of (II.3L we employ 
the powerful concept of maximal Lp-regularity (see Section 2). To this end 
we first rewrite (ED as a quasilinear evolution equation, using a smallness 
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condition like [U Assumption 2.2] which guarantees that (11.31) is a parabolic 
problem. Then we can apply the methods of the abstract parabolic theory, 

cf. ^mmm- 

Trace theory implies that the values u{t), Ut{t) at t > 0 are well-defined 
and for each t > 0 the trace map \u Ut{t))] is bounded from the class 

of solutions m into the class of the data (na. This means that the initial 
regularity is preserved for all t > 0. Concerning the asymptotic behavior of 
solutions we prove that the equilibrium u = 0 is exponentially stable. To be 
precise, we show the following result. 

Theorem 1.2. Let u be the solution in Theorenn \l.l[ Then there exist constants 
C,u! > 0 such that 

I|w(0llw2 + \\ut(t)\\^^- 2 /j, < Ce (^||uo||u /2 -I- ||wi||^ 2 - 2 /p^ , t > 0. 
Moreover, it holds that 

[t ^ e‘^^ut{t)] G W^{S, oo; Lp{n)) n W^{S, oo; W^{n)) 
for each <5 > 0 and 

^||'at(t)||w 2 + ||'att(t)||^ 2 - 2 /p) —>-0, as t ^ oo. 

This result is designed as a counterpart to m Theorem 1.3]. There, 
the authors prove that ||wtt(t)||L 2 ^ C'e“‘^‘(||Mo ||t > 0, 

which is not possible in our case, since only ui G Wp is assumed. 

However, we are able to improve the qualitative behavior of {u, Ut, Utt) with 
respect to convergence in the natural trace spaces. 

This paper is organized as follows. In Section 2, we reformulate (lESI as 
a quasilinear evolution equation of the form 

Vt + A^{v)v -I- B{v)v = F{v) in J, v(0) = vq, 

where we have introduced new functions v = [m, Ut]^, vq = [mo,wi]^. For 
sufficiently small functions z;*, we prove that the linearized problem 

Vt + A#(v.^)v + B(v^.)v = f in J, v(0) = vq, 

has the property of maximal Lp-regularity even in exponentially time-weighted 
Tp-spaces, see Theorem 12.51 

In Section 3 we prove Theorem 11.11 by applying the implicit function 
theorem in exponentially time-weighted Lp-spaces. 

Section 4 is devoted to the qualitative behavior of the solutions of dia. 
The exponential stability of the equilibrium u = 0 is a direct consequence 
of the implicit function theorem. The higher regularity of the solution u 
in Theorem 11.21 follows from the parabolicity of the problem, whereas the 
exponential decay of Utt can be seen by differentiating the equation for u and 
by studying a corresponding linear evolution equation for utt- 
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1.1. Notations 

Let p G (l,oo), J = [0,T], T > 0, and let X be a Banach space. Then 
Lp{J;X) denotes the Banach space of all p-integrable functions with values 
in X and Wp{J\X) denotes the vector-valued Sobolev space of order fc S N. 
We will frequently use the notation oWp{J;X), where 

oW^iJ-X) := {u G W^{J;X) : u(0) = 0}. 

If X and Y are Banach spaces, then C{X, Y) denotes the space of all bounded 
linear operators from X to Y. 


2. Maximal Lp-regularity for the linearized problem 

Let J = [0,r] or J = K+ := [0,oo) and let p G (l,oo). A closed linear 
operator A : D{A) —X with dense domain D{A) in a Banach space X 
is said to admit maximal Lp-regularity on J, if for each / G Lp(J;X) the 
abstract Cauchy problem 

ut + Au{t) = f{t), t G Jj u(0) = 0 

admits a unique solution u G E(J) := Wp{J]X) C\Lp{J]D{A)), cf. [2ll4l[T6]. 
Here D{A) is a Banach space for the graph norm H'lljf-l-IIA-lljf. The space 
E(J) is continuously embedded into the space trE) of bounded and 

uniformly continuous functions on J with values in the trace space 

trE = Da{1 - l/p,p) := (X, D(A))i_i/p_p, 

where (•, ■)s,p indicates real interpolation. We refer to m H [H] for more 
information on these function spaces. 

We say that the abstract inhomogeneous Cauchy problem 

Ut -b Au{t) = f{t), t G J, u(0) = uo (2.1) 

admits maximal Lp-regularity, if the solution map 

(/,uo)i->'U, Lp(J;X) X trE—> E(J) 

is a topological isomorphism. Its inverse u i—> (/, uq) will be denoted by 
{dt + A, trt=o)j where tr^^o ^ ^^(0) is the trace map. If A has maximal Lp- 

regularity on J, then also the abstract inhomogeneous Cauchy problem (12.111 
has maximal Lp-regularity on J, see [2], Section III. 1.5 and [16], Proposition 
1.2. If (12.11) has maximal Lp-regularity then the following a priori estimate is 
valid. 

II'“IIe(J) < ■^(>^) (ll^t + -b ||u(0)||trE) , uGE{J). (2.2) 

Here M{J) < oo is called a maximal Lp-regularity constant of A. It is inde¬ 
pendent of J in the following sense. For any Tq G (0, oo) there is Mq < oo 
such that (12.21) holds with M(J) = Mq for any interval J = [0,r] C [OjLq]- 
If A has maximal Lp-regularity on R_|_ we may also take Tq = oo. 
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Let us now reformulate problem (11.31) as a quasilinear evolution equation 
for the new function v = [u,ut\^. Taking into account that {u^)tt = “^UttU + 
2uf we rewrite the first equation in as 


vt ■ 


ro 

-I 1 


r 0 

01 


0 

0 

^ A 

V 

A 

0 

V = 


l — 2kvi 


L l-2kvi 


_l-2kvi_ 


where v = We now define 

A#{v) := 


m 

-I 1 

, B{v):= 

\ ® 
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\ 0 1 

0 
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AAd n 
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l — 2kvi. 




l — 2kvi 


where A u stands for the Laplace operator with homogeneous Dirichlet bound¬ 
ary conditions. This turns (USD into the quasilinear initial value problem 


Vt + A#{v)v + B{v)v = F{v) in J, v{0)=vo- (2.3) 

To ensure that the problem ()2.3D is parabolic, we concentrate on solutions v 
which are real-valued and bounded in the sense that 


sup{|ui(t,x)! : t > 0,a; G n} < l/(2fc). (2.4) 

In the remaining part of this section we investigate the linear problem 

Vt + A#{vt:)v + B{vt)v = f in J, v{0)=vo, (2.5) 

where u* is a fixed function satisfying (Ea, w : J X n —>■ is the unknown 

function, / : J x 17 —and uq : 17 —t R^ are given functions and 


A#(v^) 

A{v,,)u 


B{v,) = 


0 -I 
0 hA{vA) 

— (1 — 2k{vAji)~^ ^du. 


c'^A{v^, 


0 

0 ’ 


To shorten the notation we will sometimes drop the dependence with respect 
to u*. It is well-known that A = A(ti*) has maximal Lg-regularity in Lp(D,), 
see for instance [nisi min]. We state this as 

Theorem 2.1. Letp,q G (I,oo), fl be a bounded -domain mR", a G C{^), 
a{x) > ag >0, (cc G f7j. Then the operator Au{x) = —a{x)Au{x) with 
domain D{A) = {u € ITp(f7) : u|ao = 0} in Lp(fl), admits maximal Lq- 
regularity on M+ for each g G (1, oo). 


We conclude that A, bA, c^A have maximal Lg-regularity on R+ in 
Lp{Q) for any p,q € (1, oo). Hence also Aa := a-i- bA has this property for 
any a > 0. Consider .4^, B as operators in the Banach space Xg with domain 
Xi, defined by 

Xg := D{A) X Lp{n), Xi ■= D{A) X D{A). 

Then it follows that A^ G C{Xi,Xg) and B G C{Xg,Xg). The next goal is 
to show that 

■ 0 -I 
c^A bA • 


A I— “b ^ 
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with domain Xi in the space Xq has maximal L^-regularity on R-|.. To this 
purpose we first show in Theorem 12.21 that ^ + A has maximal L^-regularity 
on K-|. in Xq for some /r > 0. Then we use properties of the spectrum a{A) 
to prove that we may take fi = 0. 

Theorem 2.2. Let a > 0 and q G (l,oo). Then A^^a = a + A^ ■. Xi ^ Xq 
has maximal Lq-regularity on K.+ and there exists ^ > 0 such that g, -\- A = 
tL + A^ + B enjoys the same property. 

Proof. Consider the system 

vt + A#^aV = f, v{0)=0. (2.6) 

We have to show that for each / G Lq(R+; Xq) there exists a unique solution 
uGofV,^(M+;Xo)nL,(R+;Xi). 

If / = (/i) /2) and V = (ui, U 2 ), then (j2.6p is equivalent to 

dtvi + avi -V 2 = fi, 

dtV2 + AaV2 = / 2 , 

subject to the initial conditions Uj(0) = 0. Since has maximal L^-regularity 
on M+, we may solve the second equation to obtain a unique solution 

V2 €oW^iR+;Lpin))nLq(R+;DiA)). 

Inserting V 2 into the first equation we obtain the initial value problem 

dtvi+avi=g, ui(0) = 0, (2.7) 

where g ■= V 2 + fi € Lg(IR.+ ; I?(A)). It is evident that for a > 0 the op¬ 
erator {dt + a) with domain D{dt -I- a) = okTq (R-i-; I^(^)) is invertible in 
Lg(K+; D{A)) and the unique solution of (j2.7ll is given by 

^i(^) = [ t > 0. 

Jo 

This shows that A^^a has maximal Tg-regularity on R-|.. 

We finally prove the second assertion by a perturbation argument (Propo¬ 
sition 4.3 and Theorem 4.4 in 0): Since {B — a) : Xq ^ Xq is a bounded 
linear operator and a + A^ has maximal Tg-regularity on R+ , there exists 
g, > 0 such that the operator p. + A = pL + {a + A^) + {B — a) has maximal 
Tq-regularity on R+. □ 

Next, we analyze the spectrum a{A). Recall that A is defined by 

{Au){x) := —a{x)Xu{x) for u G D{A) = {u G Wp{Ll) : u\dn = 0} 

and some a G C{Cl) with a(x) > oq > 0 for all x G ^2. Since A has compact 
resolvent, the spectrum a(A) is a discrete subset of C and it consists solely 
of eigenvalues with finite multiplicity. Let A G o-(A) and u G D{A) be a 
corresponding eigenfunction. Then 

A 
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Testing the last equation with u and integrating by parts (itlan = 0) it follows 
that A is real and A > 0 by the Poincare inequality. We define 

Ai(^) := min{A : A G cr(^)} > 0, 

hence (j{A) C [Ai(A),oo). In particular, A : D{A) Lp{^l) is an isomor¬ 
phism. Suppose that A belongs to the resolvent set p{A) of A. Then the 
following resolvent formula holds. 

lx — — ~^x(A — bA) R\ 

'' ~ [I + XRx{X - bA) -XRx ■ 

Here R\ is defined by R\ := (—A^I -I- (A6 — c^)A)~^. If Xh ^ (? then 
^ ~ Xb-A {xb-A^~^) ■ 

This useful identity yields the following characterization of p{A). 

Lemma 2.3. A G piA) if and only if Xb — ^ 0 and p{X) := G p{A). 

Proof. Let A G piA). Then Rx : Lp{Q.) D{A) must be bounded. Suppose 
that Xb — = 0. Then (—A^/)“^ : Lp{Vt) —x D{A) would be bounded, but 

that is not possible. Therefore we must have Xb — (? ^0 and this implies 

Ab-c2 ^ P{^)- 

Conversely, let Xb — c^ ^ 0 and G p{A). Then Rx ; LpiVl) —x D{A) 

is bounded and we can easily check that (A—: Xq —1 is bounded. □ 


Lemma 2.4. //ReA < Aq := min{|Ai(H), then X G p{A). 


Proof. We prove the assertion by contraposition. Let a > Ai(H) and let A G C 
satisfy a = p{X) = }?/{Xb — (?). Then A has the form ^ ± ^— aA. In 

the case ^ - ac? < 0 we obtain ReA = ^ > |Ai(H). In the other case 

^ - ac > 0 we obtain that A is real and 


A > 


ab 

y 



ab^ ^ ^ 


We have shown that p{X) > Ai(AI) implies ReA > Aq. Consequently, if Re A < 
Ao, then p{X) G p{A) since cr(H) C [Ai(H),oo). By Lemma lOl this implies 
A G p{A). 

□ 


Having analyzed the spectrum for the operator A, we can continue to 
prove maximal Lp-regularity for this operator on R+. To this end, we consider 
again the linear problem (12.511 . Concerning regularity of initial data vg we 
consider the space 

:= p = D(A) X Da( 1 - l/p,p). 















S. Meyer and M. Wilke 


Note that the real interpolation space Da[^ — l/p,p) can be computed to the 
result 


Da{1 - iIp,p) 


, lip < i/2, 

{u G : u\oQ =0} , if p > 3/2, 


see e.g. [7]. In case p = 3/2 the space I?a(1 — i/p,p) looks more complicated. 

Let Y, X(J) be Banach spaces such that X(J) ^ Lijoc(</;X) and let 
a; G M. To describe exponential decay of solutions we employ the Banach 
space 


e-“X(J) := {u G Lijoc(J;X) : [t ^ e“*u(t)] G X(J)}, 


equipped with the norm ||u||e—x(j) := II [* '-t e“*u(t)]||x(j). 


Theorem 2.5. Let p G (1, oo), — Aq = sup{Re A : A G a{—A)} < 0 denote the 
spectral bound of —A and let co G [0, Aq). Then A has maximal Lp-regularity 
on ]R_|_ in the sense that 

{dt+AArt=o) : e-“(%i(R+;Xo) nLp(R+;Xi)) ^ e-“Lp(R+;Wo) x 
is an isomorphism. 


Proof. By Theorem 12.21 the operator p + A has maximal Lp-regularity on 
M+, hence A has maximal Lp-regularity on bounded intervals, which can be 
seen by multiplying the equation dtV + Av = f with By Lemma the 
spectral bound of —A is strictly negative and equals — Aq and therefore the 
spectral bound of —A+oj equals w — Aq and it is strictly negative if a; G [0, Aq). 
Then 0 Theorem 2.4] implies that A — to has maximal Lp-regularity on K-|_ 
for each w G [0, Aq), that is 

{dt + A — uj, ti’t^o) : ^p (®+; ^o) ri Lp(]R+; Xi) —> Lp(]R_|_; Xq) x Xj 
is an isomorphism. The assertion follows from [2 Proposition III.1.5.3]. □ 


3. Global wellposedness for the nonlinear problem 

In this section we prove a slightly more general form of Theorem 11.11 which 
we formulate for the transformed system (|2.3I1 : 

Theorem 3.1. Let p > max{n/2,n/4-|- 1}, w G [0,Ao) and J = K+. Then 
there exists p > 0 such that for all vg G Xj with ||uo||x., < P the problem 

vt + A{v)v = F{v), t € J, v(0) = Vo (3.1) 

admits a unique solution 

V = ivi,V2) G e-“E(J) := Xg) n Lp{J; X^)), 

such that \\vi\\buc{JxQ) < rn for some m < 1 /{2k). The solution map vg ^ v 
is a local isomorphism. 
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Theorem 11.11 corresponds to the special case w = 0. For the proof we 
will employ the following mapping property of substitution operators, which 
follows from the definition of the Frechet-derivative and the fundamental 
theorem of calculus. 


Lemma 3.2. Let k G {2, 3,...} U {oo}, T G (0, oo), J = [0, T\, and let X, Y 
be Banach spaces and let U be open in X. Suppose that S € C^(U;Y) and 
each derivative with j = 2, 3,... fc, maps bounded sets into bounded sets. 
Then the substitution operator S' := [m S' o u] belongs to the class 

C'^-\BUC{J; U); BUC{J; Y)). 


Proof of Theorem \3.1[ Step 1: In this step we prove that 


{A,F)gC°^{V;C{Xi,Xo)xXo), 
where V C X.-^ is defined as follows. Fix m G (0, (2k)~^) and let 


V := {w = {wi,W 2 ) G X^ : ||wi||oo < m} , 

then V is open in Xry, since X^ ^ C{fl) x Da{ 1 — llp,p), provided p > n/2. 

First of all, we write F in the form F{vi,V 2 ) = (0,G'i(?;i) • 02 ( 02 )), 
where Gi{vi) = (1 — 2kvi)~^, 02 ( 02 ) = 2?;| and denotes pointwise mul¬ 
tiplication. Then the maps 

Gi : {m G Loo(f^) : l|w||L,,^(n) < (2fc) t Lao(n), O 2 ■ L2p(Ll) —>■ Lp(Ll) 


are G°“-maps. Moreover, by Holder’s inequality, the pointwise multiplication 
is bilinear and continuous (and therefore 0°°) from Loo(H) x Lp(H) to Lp(Tl). 
Finally, the injection of H C Wp(Ll) x Wp~‘^^^(LI) into Loo(H) x L 2 p(Ll) is 
linear and bounded, provided that p > n/2 and p>l-|-n/4by Sobolev’s 
embedding theorem. This proves that F G 0^(V-,Xq). 

Next, we decompose 


M(z;) = Ai(v)A2 


1 0 


0 -I 

0 ^ 


bA 


V = (vi,V2) &V. (3.2) 


Here the second matrix A 2 belongs to 0°°(V-,C(Xi,Xq)) since it is inde¬ 
pendent of V and A : D(A) — Lp(Ll) is linear and bounded. The mapping 
[^1 {^1 ^ ioo(fl) : 11^111^^(0) < (2k)-^} into Loo(f2) 

and therefore [(^1,^2) is from {i;i G Loo(Ll) : ||ui||l^(o) < 

(2k)~^} X Lp(Ll) to Lp(Ll). Then the continuous embedding of Wp(Ll) into 
Lao(Ll) for p > n/2 implies Ai G G°“(H; £(Xo, Xq)). We conclude that 
MgG-(H;£(Xi,Xo)). 

Step 2 : In this step we consider the substitution operators /F : u F ou and 
A : u 1 -^ Ao u for time-dependent functions w on J x H and prove that they 
are also in appropriate spaces. Let 


i : e-‘^ (Wp(J; Xo) n Lp(J-, Xi)) e-^ BUO(J-,X^) 

be the inclusion map by [21 Theorem HI.4.10.2. & Remark HI.4.10.9.(a)]. 
Since e~‘^BUO(J; V) is open in e~^BUO(J-,X.y) it follows that 

w := {e-'^BUO(J-V)) C (Wl(J-X^) Lp(J-Xi)) = e-‘"E(J), 
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is open. We show that F G C^(W; e~‘^Lp(J; Xg)). Again we use the de¬ 
composition F(vi,V 2) = (0,Gi(vi) ■ G 2 {v 2 )) with Gi(fi) = (1 — 2kvi)~^, 
<52(^2) = 2v2 and prove that F is the composition of G°°-maps. First ob¬ 
serve that the injections 

e-‘"E( J) e-^L^{J x fl) x {e-^L^{J-, L 2 p{n)) n e-‘^Lp{J-, L 2 p{n))) 

and e~‘^LooiJ x il) ^ Loo{J x fl) are continuous and thus Then, simi¬ 
larly as in Step 1, we conclude that the map Gi : 1—>■ (1 — 2kvi)~^, 

Gi : {ui G LcaiJ X O) : ||ui||L,,^(jxn) < (2fc) L^dJ x fl) 

is G°°. Holder’s inequality implies that pointwise multiplication 
{u, u) H> MU : e~'^LodJ] L 2 p{^)) x e~‘^Lp{J-, L 2 p{fl)) -?■ e~‘^Lp{J; Lp{i})) 
is continuous. We conclude that 

G2 : 1^2 ^ 2vl e-‘^L^{J-,L 2 p{n))ne-‘^Lp{J-L 2 p{n)) -G e-^Lp{J x H) 
is G°°. Since pointwise multiplication is also continuous as a mapping 
Loo{J X H) X e~^Lp{J X H) —>• e~‘^Lp{J x H), 

we obtain that F -.W ^ e~^Lp{J\Xg) is G°°. 

In Step 1 we have shown that A G G°°{V] C{Xg,Xi)). The previously 
used decomposition of A implies that all its derivatives are bounded on V. 
From Lemma and the embedding e~^BUG{J\Xd ^ BUG{J] Xj) we 
infer that 

AgG^ {e-‘^BUG{J; V);BUG{J;C{Xi,Xo))) . 

Here we use the same notation for A and the substitution operator induced 
by A. Furthermore it is easy to see that the mapping 

[{B,v) ^ Bv] : BUG{J]C{Xi,Xg)) x e"“Lp(J; ATi) e“‘^Lp(J;Xo) 

is continuous and bilinear, hence G°°. Defining G{v) := A{v)v it follows that 
GGG^{W-e-^Lp{J-Xg)). 

Step 3: Define a mapping H :W x X^ —> e~^ Lp{J\ Xg) x X^ by 

H{v, Vg) := {Vt + A{v)v - F{v), tlt=g V - Vg). 

By Step 2, we may conclude that 

H G G\W X XY,e-‘^Lp{J;Xg) x X^, 

since dt and tr^^o are linear operators on W. 

Note that u € IF is a solution of (EH) if and only if F[{v,vg) = 0. 
Obviously it holds that iL(0,0) = 0 and 

DyH{Q, 0)ru = {wt + A(0)m;, tr^^o w), w GW. 

By Theorem 12.51 the operator 

(5t + A(0),trt=o) : e-^ Xg) ^ Lp{J-Xd) ^ e-^Lp{J-Xg) x X^ 

is an isomorphism with bounded inverse, hence by the implicit function the¬ 
orem there exists an open ball ]Bp(0) C X.y and a unique function </> G 
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C'^(]Bp(0); such that H{<^{vo),vo) = 0 for all vq G Bp(0) and (/)(0) = 
0. Therefore v := (j){vo) G IT is the unique solution of (13.11) and i{v) G 
e~‘^BUC{J]V), hence 111^1(^)1100 < m < l/(2fc) for all t G J. □ 

4. Exponential Stability 

Let 

E(M+) :=WliR+;Xo)nLp{R+;Xi). 

By Theorem 13.II with T = oo there exists p > 0 and a function (j) in the class 
C'^(Bp(0); e“‘^E(]R+)) with ^(0) = 0 such that v = cj){vQ) G e“‘^E(]R_|_) is the 
unique solution of (EH). Then we obtain 

4)iyo) = j ^^0, 

for all vq G ]Bp(0). Since (j) is C^, we see that for each e > 0 there exists 
Po = Po(0) > 0 such that 

||</''(wo) - <(''(0)||£(Xo;e—E(R+)) < £, 

whenever ||uo||x.^ < Po < P- It follows that 4)'(vo) is uniformly bounded in 
£(X.y; e““E(R+)) w.r.t. vq if ||wo||x.^ < Po < P- Therefore we obtain the 
estimate 

e—'E(R+) < C'lkollx., 

which is valid for alH > 0, |fo|xo < Po < P and C > 0 does not depend on 
Vo and t. This implies exponential stability of the equilibrium 0 in Xj and 
the first estimate in Theorem o If w = 0 then we still have ||u(t)||xo —t 0 
as t —>■ oo, since v G E(IR+) ^ C'o(K+;^7). 

But we can show even more. By [H Theorem 5.1] we obtain that 

veWp^(S,T;Xo)nWp\S,T;Xi), (4.1) 

for all 0 < 5 < T < oo, since (A,F) G C^(V; £(Xi, Xq) x Xq) and since 
A(v(t)) has maximal Lp-regularity for each fixed t G [0,T], by Theorem 12.51 
(recall that ||£i(t)||oo < m < l/(2fc) for all t > 0). In particular, (14.11) implies 
that Vt{t) G X^ for each t > 0. 

We will show that we may even set T = oo in EH and that e‘^*Vt{t) —>• 0 
in Xj as t oo. The result on higher parabolic regularity enables us to 
differentiate equation dUD w.r.t. t > 6 to the result 

vtt(t) + [A'{v{t))vt(t)]v(t) -b A{v(t))vtit) = F'{v(t))vtit), t > S. 

We now distinguish between the function w \= Vt and the fixed solution v, 
hence fo is a solution to the linear problem 

wt + A{t))w{t) = B{t)w{t), t > 6, w{S)=vt{S), (4.2) 

where 

B(t)w(t) := F'{v{t))w{t) + [^(0) — A{v{t))]'w{t) — [A'{v(t))w(t)]v(t). 
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The operator B{t) is well-defined, since v{t) £ Xi = D{A) x D{A) for a.a. 
t > 0. 

We claim that for each e > 0 there exists a sufficiently large 5 > 0 such 
that the estimate 


||i3(-)w;||e— Lp( 5 ,oo;Xo) < II e-"E( 5 .oo) (4-3) 

holds for each w £ e“‘^o]E((5, oo), where olE(d, oo) := {u £ E((5, oo) : u{5) = 0}. 
Since e^*v{t) —)• 0 in as t ^ oo and v{t) £ V for all t > 0, the desired 
estimate follows for the first and the second term, so we concentrate on the 
third one. We use that the spaces e“‘^E(d, oo) are continuously embedded into 
e~‘^BUC{6,oo-,Xy) and the embedding constant is independent of 6, which 
follows from identities like 


ll^^lle-“E(5.oo) = l|[i '->■ + (5)]||e(0,oo)- 

From H3.2D we obtain 


\\[A'{v)w]v\\e-<-L^(^S,oo-,Xo) 


2k 


+ iAv,) 


e~^ Lp{5,C!0‘,Lp{Q.)) 


Since ||ui(t)||oo < m < l/(2fc) for all t > 0 and w £ e ‘^o]E((5, oo) ^ 
e~‘^BUC{6,oo-,X-y) and D{A) ^ Loo(fl) it follows that 

||[^'(u)w;]?;||e-^i^(5_oo;Xo) 

< C'm||w||e-‘^B;7C(5.oo;X.^)||(c^^'yi + ^^^ 2 ) || e-"Lp(5,oo;Lp(n)) 

^ Cm ||^||e“‘^E({5,oo) II II e“^Lp ({5,oo;Xi) • 

Since v £ e““E(IR+) is the solution of (EU and thereby a fixed function, 
it follows that ||r’||e-‘^Lp((5,oo;Xi) < ^/Cm, provided that d > 0 is sufficiently 
large. This yields the claim. 

Finally we write 

dt + M(0) - B{t) = [I- B{t){dt -t- ^(0))-i](at + M(0)). 


We show that the operator on the left side is an isomorphism from e“‘^oE(^! c») 
to e~‘^Lp{S,oo;Xo). To this end, it suffices to show that the first oper¬ 
ator on the right hand side is an isomorphism from e~‘^Lp{S, oo] Xq) to 
e~‘^Lp{S,oo;Xo). We have 

||5(-)(9t -I-^(0)) < e||(9t -f M(0)) ^/l||e-‘^E((5,oo) 

< eC'||/l||e-^i^(5_oo;Xo)- 

Here the maximal regularity constant C > 0 of (9t + M(0)) does not depend 
on (5 £ R+. To see this, let denote the norm of {dt — uj + ^(0))“^ : 
Tp(R+;Xo) —>■ E(R+). Using a translation from (<5, oo) to (0,oo) we obtain 

l|■u||e-“E(^.oo) = Il[i + (5 )]||e(r+) 

< M^||[< (9t - w + A{0)){e‘^^*^^'’u{t + d))]||Lp(R^;Xo) 

= Mi^\\{dt + -4(0))u||e-u.j^p(,5_oo.Xo). 







Regularity and long-time behavior for the Westervelt equation 13 


If we choose eC < 1, hence 5 > 0 sufficiently large, a Neumann series argu¬ 
ment implies that the operator dt + A{0) — B{-) is invertible from e“‘^o]E((5, oo) 
to e~‘^Lp(6,oo; Xo) with bounded inverse. 

From now on we fix such a 5 > 0, say (5* > 0, and define w{S*) := 
Vt{S*) G X^. By Theorem l2.5l the function w^(t) := belongs 

to e“‘^E(^*, oo). Then the shifted function w := w — = Vt — 'w^, solves the 

problem 

wt{t) + A{0)w{t) — B{t)w{t) = t>S*, w(S*) = 0. (4.4) 

It is easily seen that B{-)w^, is in e~‘^Lp{S*,oo-, Xq). By the preceding argu¬ 
ments it follows that w G e“‘^o]E((y*, oo) or equivalently vt G e“‘^E((5*, oo). 
Property (ED with T = S* thus yields Vt € e ‘^E((5, oo) for each S > 0 (the 
exponential weight plays no role on the bounded interval (d, 5*)). This implies 
in particular that e‘^*Vt(t) — > 0 in X^ = D{A) x Da^A ~ 1/PiP) as t —t oo. 
We summarize these results in 

Theorem 4.1. Let n G N, p > max{n/2, n/4-|- I}, lo G [0,Ao) and let vq 
satisfy the conditions of Theorem 1,9. IL Then the equilibrium v^, = 0 of (EU) 
is stable in X^ and e“*||?;(t)||x.^ —> 0 as t —> oo. Moreover, it holds that 
Vt G e“‘^E(5, oo) for each S > 0, hence e“*||ut(t)||x., —t 0 as t ^ oo. 
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